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Abstract 

Dynamical properties are studied for escaping particles, injected through a hole in an oval billiard. The dynamics is considered for 
both static and periodically moving boundaries. For the static boundary, two different decays for the recurrence time distribution 
were observed after exponential decay for short times: A changeover to: (i) power law or; (ii) stretched exponential. Both slower 
decays are due to sticky orbits trapped near KAM islands, with the stretched exponential apparently associated with a single group 
of large islands. For time dependent case, survival probabiUty leads to the conclusion that sticky orbits are less evident compared 
with the static case. 
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1. Introduction 

Billiards are dynamical systems where an ensemble of mov- 
ing particles do not interact with each other and suffer specular 
ireflections with the boundary y] |2|] . Applications of billiards 
can be made to different physical systems including experi- 
ments in reflection of light from mirrors [3], superconducting 
i^ and confinement of electrons in semiconductors by electric 
potentials [5, o], wave guides |7], microwave billiards f8l|9[], 
ultra-cold atoms trapped in a laser potential [10, 11, 12, 13|1 and 
ialso mesoscopic quantum dots 01411 . 

' Billiard dynamics falls into three main classes namely: (i) 
jregular ||2|; (ii) intermittent |15| or; (iii) totally chaotic behav- 
ior 01611 . Two examples of case (i) are the circular billiard [Uj], 
which preserves both the energy and angular momentum, and 
the elliptic billiard which preserves the energy and the product 
pf the angular momenta about the two foci yfl . Case (ii) holds 
for typical billiard shapes, including many models considered 
jso far il7i, 18, 19, 201 , producing a mixed phase space struc- 
ture in the sense that elliptic islands with fractal boundaries, 
generally surrounded by a chaotic sea that often is confined by 
invariant spanning curves, can all be observed. The latter case 
(iii) includes Sinai 112 ill dispersing billiards with comers as the 
diamond [22], or periodically extended such as the finite hori- 
zon Lorentz gas [23, 24]. Billiards may also be constructed 
that are fully chaotic according to the defocusing mechanism 
lll6ll25n . However many famous examples such as the Sinai 
billiard consisting of a square with a circular obstacle (or the 
equivalent infinite horizon Lorentz gas 11261 12711 and the Buni- 
movich stadium have a regular family of periodic orbits, thus 
making the dynamics intermittent as in case (ii), even though 
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the phase space is a single chaotic ergodic component. This 
makes them analytically tractable models for mixed systems, 
similar to the more recently introduced mushroom billiards, 
which have chaotic and regular regions of phase space, but sep- 
arated by smooth (rather than fractal) boundaries |J28| |29|. 

A time perturbation to the boundary may be considered, for 
example due to thermal vibrations in solids [30], with amplitude 
and typical frequency related to the temperature. Depending on 
the type and shape of the billiard, such a time-dependence leads 
to the so called Fermi acceleration (FA) [31]. This phenomenon 
consists in the unlimited energy growth of the bouncing particle 
upon collisions with the, presumably, infinitely heavy moving 
boundary. Several diff'erent kinds of perturbation can be con- 
sidered in different billiard-like models ll32[ l33[ l34l l35[ 13611 . 
As claimed in the Loskutov-Ryabov-Akinshin (LRA) conjec- 
ture 113711 . if the dynamics of the particle is chaotic while the 
boundary is static, thus this is a sufficient condition to observe 
FA when a time perturbation to the boundary is introduced. Re- 
cently it was shown 1 38] that even a time-dependent elliptic bil- 
liard, which is integrable for the static boundary, can also gener- 
ate FA thanks to the appearance of a stochastic layer replacing 
the separatrix curve in the phase space. Moreover, the existence 
a heteroclinic orbit could extend the LRA conjecture [39] and 
the unlimited energy growth can be observed even in (some) 
integrable billiards. The occurrence of an exponential FA was 
reported in a time varying rectangular billiard [l40ll . which was 
latter explained [41] as due to a sequence of highly correlated 
motion which consists of alternating phases with free propaga- 
tion motion along the invariant spanning curves of the Fermi- 
Ulam model; see Ref. [42[ for the localization of such curves 
for a family of mappings whose angle is a diverging function of 
the action in the limit of vanishing action, including the Fermi- 
Ulam model. 

When a hole is introduced in the boundary therefore letting 
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the particle leave the billiard region, we may consider related 
problems of escape (initial conditions in the billiard) or recur- 
rence (initial conditions at the hole); here we consider recur- 
rence, but virtually all of the following discussion applies also 
to escape, with a different exponent in the distributions if there 
is a power law decay. It is known that, for fully chaotic dy- 
namics, the recurrence time distribution, i.e. the time the par- 
ticle stays confined in the billiard domain, is characterized by 
an exponential decay |43]. On the other hand, for intermittent 
including mixed phase space dynamics where there is sticki- 
ness, generated from a finite time (but arbitrarily long) trapping 
near periodic/elliptic regions, a power law decay is observed 
1 4411 . Recently, an investigation of a mushroom billiard led to 
the characterization of families of marginally unstable periodic 
orbits [29] responsible for trapping the particle in sticky do- 
mains, including their effects on the escape problem. More- 
over, for a time dependent potential well ll45ll . the dynamics of 
the particle is shown to be fully chaotic for the low energy do- 
main and reaching elUptic islands as far as the energy increases 
until finding a limitation marked by the existence of an invariant 
spanning curve. For the time dependent potential well, a hole 
in the energy space was introduced letting the particle escape. 
Therefore for the low energy regime an exponential decay was 
observed while a slower decay characterized either as a power 
law or stretched exponential marks the regime of higher energy 
and consequently long time. Thus opening a billiard by consid- 
ering particles escaping through a hole is a good means of iden- 
tifying and describing various kinds of intermittency present in 
the dynamics. 

The oval considered here is defined by a finite Fourier series 
in polar coordinates, and the mixed phase space of oval billiards 
was first described by Berry in 1981 lllSll . Since then it has re- 
mained a popular example of a billiard with mixed phase space, 
including for generalizations to time dependent boundaries [46] 
and wave chaos in theory ll47ll and microresonator experiments 
1 4811 . In this paper we revisit the oval billiard considering both 
the static and time-dependent boundary. We consider a hole for 
the first time, seeking to understand and describe some proper- 
ties of particles returning to the hole. 

This paper is organized as follows. In Sec. |2] we describe 
the model with fixed boundary, detailing results for the recur- 
rence times. Section [3] considers the moving boundary where 
the equations of the mapping are derived. The results for the 
recurrence times are discussed here also. Final remarks and 
conclusions are presented in Sec. |4] 




Figure 1 : Illustration of the angles that describe the dynamics and an escaping 
trajectory. 



ter e = the circular billiard is obtained leading to a foliated 
phase space [IJ. Therefore chaos is not observed. In the case 
of e ^ but considering e < e^ = l/(p^ + 1), the billiard is 
convex, and the phase space contains both elliptic islands, in- 
variant spanning curves corresponding to rotating orbits (also 
called whispering gallery orbits) and chaotic regions [49 1 while 
for e > Ef the billiard is no longer convex; all the invariant tori 
are destroyed 1,19.1 however some elliptic islands survive. 

The dynamics of the particle is described by a two- 
dimensional nonlinear area preserving map T for the vari- 
ables (9„,a„) where 6„ is the angular position of the parti- 
cle and a„ is the angle that the trajectory of the particle does 
with respect to the tangent vector of the boundary at the an- 
gular position 9n (see Fig. [T]i. The index n corresponds to 
the n''' collision of the particle with the boundary. Using po- 
lar coordinates one has that X{9n) - \l + ecos(/90„)] cos(6'„) 
and Y(9n) - [1 -H ecos(/:i0„)] sin(6'„). For an initial condition 
(9„,a„), the angle between the tangent and the boundary at 
the position X{9„) and Y(9„) with respect to the horizontal is 
0„ - aictan[Y'(9n)/X'(9„)]. Between collisions, the parti- 
cle travels with a constant velocity along a straight line until 
reaches the boundary. The equation that gives the trajectory of 
the particle is 



F(0„+i) - Y(9„) = tan(a„ + cf>„)[X(9„^i) - X(9„)l 



(2) 



2. The static oval billiard, the mapping and escaping parti- 
cles 

The model we consider in this section consists of a classical 
particle confined to move in a domain which the radius of the 
boundary is given by the following equation in polar coordi- 
nates 



R(9,e,p)^ 1 +ecos(p9) 



(1) 



where e is the amplitude of the circle's perturbation, 9 is the 
angular coordinate and p > is an integer For the parame- 



where (p„ is the slope of the tangent vector measured with re- 
spect to the positive X-axis, X(9„+i) and Y{9„+i) are the new 
rectangular coordinates of the collision point at 9„+i, which is 
numerically obtained as solution of Eq. (|2]i. The angle between 
the trajectory of the particle and the tangent vector to the bound- 
ary at 9n+i is 



««+! = 0«+l - (a„ + 4>n)- 



(3) 



Figure \T\ illustrates the corresponding angles and a escaping 
particle from the billiard. The mapping that describes the dy- 
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Figure 2: Phase space for the oval billiard for the control parameters: (a) e = 
0.07 and (b) e = 0. 1 . (c) shows a typical periodic orbit and (d) a zoom-in of a 
specific region of (b). 



namics of the model is thus given by 



( F(0„+i) = /?(0„+i)sin(0„+i) - Yi0„)- 
T : I tan(a„ + 0„)[^(0„+i) cos(0„+i) - X(6„)], 

[ a„+i - (pn+i - (an + 4>n) 



(4) 



where 0„+i is numerically obtained as solution of 
F(0n+i) = with R(0„+i) = 1 + ecos(p0n+i) and 
4>„+i =arctan[r(0„+i)/r(0„+i)]. 

A typical phase space for the static version for different con- 
trol parameters together with a visualization of a period three 
orbit is shown in Fig. |2] The parameters used in the figure were 
p = 3 and: (a) e = 0.07 < e,., (b) e = 0.1 = e^. Figure He) 
shows a period three orbit indicating corresponding region in 
the phase space of (a) while (d) shows zoom-in of a region near 
a elliptic island of (b). 

Let us now consider that the boundary has a hole through 
which the particles are injected and can escape, as shown in 
Fig. [U We assume the hole is localized in 6 (0, h) where h 
is a parameter We simulated different values of h < tt/IO how- 
ever in this paper we fix it at /i = 0.1. The procedure used to 
consider the escape of the particles assumes the evolution of an 
ensemble of initial conditions. Indeed we consider 10'' differ- 



Figure 3: (a) Histogram of frequency for the escaping orbits from the billiard, 
(b) Corresponding survival probability, obtained by integration of the histogram 
shown in (a). The parameters used were p = 3 and e = 0.07 < er> f = 0. 1 = e^ 
and e = 0.13 > fc- 



ent initial conditions in a window where 10^ 00 are uniformly 
distributed along 0o e {0,h) while a window of 10^ different 
Q-Q also uniformly distributed along cq e (0, n). Each one of 
them was let to evolve a maximum of 10^ collisions with the 
boundary, if it did not escape before. When the particle reaches 
the region of the hole for the first time, the number of collisions 
with the boundary spent up to that point is registered and the 
particle is assumed to escape. A new initial condition is then 
started and the procedure is repeated until all the ensemble is 
exhausted. The histogram of frequency of escaping particles, 
represented as H{n) is shown in Fig. [3la) for three different pa- 
rameters, as labeled in the figure. For a fixed p = 3, the param- 
eter e = 0.07 < Ee causes the phase space to have both elliptic 
islands as well as invariant spanning curves lll9ll correspond- 
ing to the so called whispering gallery orbits. The presence 
of the elliptic islands leads the dynamics of some initial con- 
ditions to experience a sticky behavior that can be long. The 
invariant spanning curves are destroyed for the cases of e = 0.1 
and e - 0.13. Figure [3j a) shows the behavior of the histogram 
of escaping orbits. The horizontal axis denotes the number of 
collisions the particle suffered with the boundary before escap- 



ing while the vertical one corresponds to the fraction of orbits 
which escaped at the n* collision. Given that p -3 one can see 
that the histogram shows a lower value for three bounces with 
boundary as compared with n -2 and n - A. This reduction is 
related to the stability of period three orbits in the phase space, 
therefore trapping the particle close to this region. For large 
n we see a long tail which corresponds to sticky orbits. The 
integration of such histogram gives the so called cumulative re- 
currence time distribution, which is defined as 



1 ^ 



A^rec(n) : 



(5) 



where the summation is taken along the ensemble of N - 10^ 
dififerent initial conditions. A^rec(n) denotes the number of initial 
conditions that do not escape through the hole (ie recur) until a 
collision n. When Eq. Q is evaluated in a fully chaotic dynam- 
ics its behavior is an exponential [43] while for a mixed phase 
space where intermittent orbits exist along the phase space a 
power law is observed 04411 . We have shown recently 115 Oil that 
the existence of elliptic islands may also lead to a stretched ex- 
ponential decay. Figure|3lb) shows the behavior of three curves 
of P(n) vs n for the same set of control parameters used in 
Fig. |3ja). For e - 0.07 < e^ the decay is exponentially fast 
at the beginning until about 200 collisions of the particle with 
the boundary when the curve changes to a slower decaying 
regime marked by a power lay with exponent -2.431(4). For 
e = 0.1 - ec, the invariant spanning curves creating the whis- 
pering gallery orbits are destroyed. The decay of P at the be- 
ginning is the same for e - 0.07 when a hump appeared around 
n = 500 lasting until n = 1500. Indeed the hump is described 
by a stretched exponential of the type 



P = Po exp(/7 n'^) 



(6) 



with the coefficients Pq - 0.0054, b - -0.104 and y - 
0.4674 = 0.5. From the 10^ different initial conditions, the re- 
gion corresponding to the hump is due to 947 initial conditions. 
The major part of the trapping happens near a period three orbit 
as shown in Fig. 15 a) with the corresponding sticky orbit plot- 
ted in Fig. SJb). Finally, for e = 0.13 > fc the elliptic regions 
in the phase space are reduced and the recurrence distribution 
decays rapidly at first. After n > 1360, some of the initial con- 
ditions are trapped for long time in a sticky region. The few 
orbits trapped for a long time were mostly observed near a pe- 
riod twelve orbit, as shown in Fig. |4jc) with its corresponding 
plot in the phase space shown in Fig. |4jd). 

3. Time dependent oval billiard and escaping particles re- 
sults 

This section is devoted to discussing the recurrence time dis- 
tribution of the time-dependent oval billiard. We first construct 
the mapping that gives the precise description of the dynamics. 
The radius of the boundary in polar coordinates, to include the 
time-dependence, is now written as 



(a) 



(c) 




Figure 4: Plot of a typical long-lived orbit in the billiard (a,c) and its cor- 
responding representation in the phase space (b,d) for the parameters p = 3 
h = 0.1 and: (a,b) f = 0.1 = e^; (c,d) e = 0.13 > e^. 



where a corresponds to the amplitude of oscillation of the 
boundary. The introduction of the time perturbation to the 
boundary produces two new additional variables that have to 
be considered, namely: (i) the velocity of the particle, V, and 
the time t. The map describing the dynamics has now four dy- 
namical variables, i.e. T{6„,a„,V„,t„) = (0„+i,a„+i,V„+i,tn+i). 
Supposing the initial conditions (6„, a„, V,,, f„) are given, a sim- 
ilar procedure as made in the previous section can be used to 
describe the position and trajectory of the particle. Then the 
instant of the collision is obtained by the numerical solution of 
the following equation 



Rn(0, t) = RM t) 



(8) 



Rb(6, f) = 1 -H e[l + fl cos(r)] cos(p6) 



(7) 



where Rh(0,t) - 1 + e[l + flcos(f„ + t)]cos{pOp) and 

Rp(t) - Jxj(t) + Yj(t) with the corresponding angle 0p = 

arctan[Yp(t)/Xp(t)l Xp(t) = X(6„, t„) + \V„\ cos(<f>„ + a„)(t - t„) , 
and Ypif) = Y{9„, t„) + |y„| sin(0„ -i- a„){t - t„) , with t > t„. 

The angular coordinate at the new collision, 0„+], is numeri- 
cally obtained from Eq. dHJ via a numerical procedure similar to 
the molecular dynamics method leading to an accuracy of 10 '^ 
in the time of the collision. Given 6„+i, the time at the collision 




Figure 5: Plot of five snapshots of a trajectory for an escaping particle in the 
time varying bilUard. The control parameters used, for visual | 
p = 3, 6 = 0.1 and a = 0.8. 



I purposes, were 



IS written as 



tn+l — f n + 



[(AXy + (AYY] 
v„ 



211/2 



(9) 



where AX = Xp(6„+i)-Xp(e„) and AF = Yi,(e„+i)- Yp(6„). The 
velocity of the moving boundary is given by 



Vb(tn+i) = -easin(f„+i)cos(/?6l„+i) 
X [cos(6'„+i)f + sin(6'„+i)/] 
The reflection laws used are 



(10) 

(11) 
(12) 



where the upper prime denotes the variables are represented in 
the moving referential frame. From Eq. (ITTt one concludes 
that the tangent component of the velocity does not indeed suf- 
fers any modification after the impact. Returning to the inertial 
frame of reference, we obtain that 

v„+i • T„+i = v„[cos(q'„ + ^„) cos{(p„+i)] + 

+ v„[sin(Q'„+0„)sin((;A„+i)] . (13) 

Considering Eq. (fT2l i. in the rest referential frame, the normal 
component of the velocity of the particle is 

N„+i = -v„[-cos(q;„ + <;A„)sin(0„+i)] + 



v«+i 



+ v„[sin(Q'„ + (p„) cos(<^„+i)] - 

- lea sin(f„+i) cos(p0„+i)[- cos(6l„+i) sin(0„+i)] 

- 2eflsin(f„+i)sin(6'„+i)cos(0„+i)] . 



The velocity of the particle immediately after the impact is 



v„+i = -\/(^Wi ■ r„+i)2 + (v„+i ■ N„+i)^ 



(14) 



and finally, the angle that the particle leaves the boundary, mea- 
sured with respect to a tangent to the point d„+i is written as 



a„+i = arctan 



Vn+l ■ N„- 



Vn+l ■ T„+i 



(15) 
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Figure 6: (a) Histogram of frequency for the escaping orbits from the billiard, 
(b) Corresponding survival probability, obtained by integration of the histogram 
shown in (a). The parameters used were p = 3 and e = 0.07 < er> f = 0. 1 = fc 
and € = 0.13 > e^ and a = 0.1. 



For fl ?i the particle can gain or lose energy upon colUsions 
with the boundary and given that the phase space has chaotic 
components, unlimited energy growth is observed plT there- 
fore confirming the LRA conjecture. Figure |5] shows 5 snap- 
shots of an orbit as well as the corresponding position of the 
wall at the instant of the collisions. The parameters used, only 
for visual purposes were p - 3, e -0.1, a = 0.8 with /i = 0.1. 

The histogram of frequency for the escaping orbits is shown 
in Fig. |6la) for the parameters /:> = 3, a = 0.1 and the same 
three different e, as used in the static case. We can see from the 
histogram (Fig. |6la)) that the escaping particles at 3 collisions 
are still less observed than the ones for 2 and 4 collisions. The 
survival probability was considered for different values of e and 
for three different values of initial velocity. For e = 0.07, we 
can see in Fig. |6jb) that the survival probability decays expo- 
nentially fast and few orbits keep trapped at large n therefore 
leading the curve to slower the decay at the end. The slower 
initial velocity Vq = 0.1 seems to affect less the dynamics while 
a short tail of slower decay is observed for Vo = 1 and Vq = 10. 
For e = 0.1 (see Fig. [SJc)), the survival probability for Vo = 10 
is marked by a hump starting at n ^ 400 while few orbits are 



trapped in sticky domain for Vo = 1 and Vq = 0.1. On the other 
hand, when e = 0.13 (see Fig. |6ld)), no significant changes 
from fast exponential decay was observed as dependent on the 
initial velocity. 

Let us discuss the results obtained for the survival probabil- 
ity considering the two cases of (i) static and (ii) time-varying 
boundary. For the static boundary, the velocity of the particle 
is constant and for a mixed phase space where invariant span- 
ning curves and KAM islands coexist, the sticky behavior ob- 
served is larger as compared to the case where invariant span- 
ning curves are destroyed (e > e^ (see Ref. |119|])) and chaotic 
sea is limited by KAM islands only. Therefore the decay of 
the survival probability starts exponentially at short collisions 
and suddenly it is marked by a changeover turning into a power 
law for large number of collisions. On the other hand, when 
the invariant spanning curves are destroyed but the period three 
region in the phase space still influences the dynamics, several 
instances of trapping were observed leading the dynamics to 
spend long time near period three orbits. The decay of the sur- 
vival probability starts exponentially fast at short coflisions and 
suddenly it changes to a slower decay being characterized by a 
stretched exponential. The decay is slower than exponential but 
is still faster than a power law. Indeed, as the control parame- 
ters are varied and invariant spanning curves are destroyed, one 
can observe a continuum spectrum of decay ranging from expo- 
nential to a power law. This variation is still an open problem 
and extensive theoretical and numerical investigations should 
be made to describe it properly. As the parameter e rises, the el- 
liptic region in the phase space decreases and the exponentially 
fast decay of the survival probability is most evident. However 
trapping is still observed for large times. In our case we ob- 
served a few long orbits trapped near a region of period twelve. 
Such orbits indeed slow the decay of the survival probability at 
the very long time but were observed only in a few trajectories. 

For case (ii) where a time varying boundary is considered, 
the velocity of the particle is no longer constant. The LRA 
conjecture claims [37], the chaotic dynamics of the particle for 
the static boundary leads to the unlimited energy growth when a 
time perturbation is considered; numerical studies of this model 
are consistent with this prediction [46, 51]. A consequence is 
that a particle with high energy collides many more times with 
the boundary in a given interval of time while compared with a 
lower energy particle at the same interval of time. Over a small 
number of coflisions the billiard it sees is effectively static, and 
it is likely to escape well before Fermi acceleration is evident. 
This is observed particularly for the case of e = 0. 1 and Vq = 10 
where a hump in the survival probability is evident. 

4. Concluding remarks 

We have studied some dynamical properties of an oval-like 
billiard with a hole in the boundary, considering both static as 
well as time dependent boundaries. For the static case, the 
recurrence time distribution of the hole has a fast decay for 
short collisions changing the decay either to a power law or 
stretched exponential, depending on the control parameter. The 
power law observed for e = 0.07 in the static case has slope 



-2.431(4) while the stretched exponential for e = 0.1 is given 
by P = Po exp(fe «'') with coefficients Pq = 0.0054, h = -0.104 
and J - 0.47. The sticky orbits present in the dynamics are 
responsible for slowing the decay of the recurrence time dis- 
tribution. For the time dependent case, the survival probability 
was considered for different values of e and for three different 
values of initial velocity. For e = 0.07, the lower initial veloc- 
ity seems to affect less the trapping orbits while a short tail is 
observed for larger initial velocities indicating a sticky regime. 
For e = 0.1, the initial Vq = 10 lead the survival probability to 
exhibit a hump starting at n x; 400 therefore indicating a sticky 
regime. Indeed at that large energy, the particle suffers many 
more collisions with the boundary at the same interval of time 
as compared to a low energy particle, hence seeing less the in- 
fluence of the moving boundary compared with a lower energy 
particle, therefore seem to be more susceptible to sticky behav- 
ior. For e = 0.13 the initial velocities considered do not seem to 
change significantly the fast exponential decay as observed for 
e = 0.1 and e = 0.07. 

The observation of stretched exponential decays, as in 
Ref. ]'50J, invites further investigation, as previous studies have 
concentrated on algebraic decay models. We would like to re- 
mark that as in Ref. ll50ll . a stretched exponential decay is ob- 
served where there is a single prominent elliptic periodic orbit 
responsible for the stickiness, and also that the exponent j is 
very close to 1/2. In general, opening a dynamical system with 
a hole is a very effective method of elucidating the structure of 
intermittent dynamics. 
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